Non-Ohmic variable-range hopping transport in one-dimensional conductors 
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We investigate theoretically the effect of a finite electric field on the resistivity of a disordered one- 
dimensional system in the variable-range hopping regime. We find that at low fields the transport is inhibited 
by rare fluctuations in the random distribution of localized states that create high-resistance "breaks" in the hop- 
ping network. As the field increases, the breaks become less resistive. In strong fields the breaks are overrun 
and the electron distribution function is driven far from equilibrum. The logarithm of the resistance initially 
shows a simple exponential drop with the field, followed by a logarithmic dependence, and finally, by an inverse 
square-root law. 

PACS numbers: 72.20.Ee, 72.20.Ht, 73.63.Nm 
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Motivation and main results. — Potential applications in nan- 
otechnology has generated a widespread interest in one- 
dimensional (ID) conductors, including carbon anotubes, 
nanowires, and conducting molecules. Before such applica- 
tions may become practical, a number of fundamental physics 
questions need to be resolved. One of them is the effect of 
disorder on transport properties. Despite a great progress in 
nanofabrication, the disorder is hard to get rid off, especially 
if a low-cost processing is desired. The disorder causes local- 
ization of electron states, so that the low-temperature transport 
in ID systems is often of the hopping type |2D - That is why 
a new look at the theory of a ID hopping conduction may be 
timely. One of the open questions here is the mechanism of 
the non-Ohmic hopping transport, observed in a number of 
recent experiments 1^,0,0]. Especially interesting and chal- 
lenging is the ultimate low-T limit where the variable-range 
hopping (VRH), i.e., hopping between distant localized states 
dominates. Remarkably, this problem has long remained un- 
touched by theorists except for one early numerical study jjj. 
Recently, however, the following analytical formula has been 
proposed l ioll : 



p oc exp(Ck B To/Fa) 
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where p is the resistivity, F — e'E is the product of the elec- 
tron's charge and the electric field, a is the localization length, 
To a parameter of dimension of temperature defined below, 
and C is a numerical coefficient left undetermined. This for- 
mula has immediately become a matter of controversy be- 
cause in its derivation a very important and nontrivial peculiar- 
ity of the ID geometry was overlooked: unavoidable highly 
resistive segments ("breaks") on the conducting path j^ElE!- 
In an attempt to include this physics, Eq. Q was subsequenty 
revised in a follow-up article [ 1 1 J . Yet the actual situation 
is more intriguing. As demonstrated below, the problem has 
quite a rich regime diagram. The resistivity may take three 
different functional forms, depending on the strength of the 
electric field (weak, intermediate, or strong). It is Eq. Q that 
proves to be correct in strong fields, F 3> ksT /a, not its sub- 
sequent revision ITUl . We compute the exact coefficient in 
Eq. <□ to be C = 8. The main message of our analysis is that 
the role of the breaks diminishes and eventually becomes in- 



significant as the electric field increases. We hope that these 
ideas and new results may provide a much needed theoretical 
basis for the purposeful experimental study of the non-Ohmic 
VRH expected to be reported in the near future. 
Model. — Since the problem has proved to be rather delicate, 
in this study we approach it within a simplified model of non- 
interacting electrons. (In real systems interaction effects can 
be suppressed by placing the ID conductor near a metallic 
gate.) The localized states are regarded as random points 
{xj,Ej) with the Poisson distribution of a constant density g 
in the position-energy space. This is justified at low T where 
electron transport is determined by a narrow range of ener- 
gies. Similarly, at low enough T we can neglect the field- 
dependence of the localization length a because the electric 
fields of interest, say, F ~ ksT/a, are small. 

The rate of hopping F n from site j to site k is given by the 
usual model form 111 11311 



=fj(l -/,)exp 
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max(D jk ,0) 
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where y of dimension of frequency is determined by the 
electron-phonon coupling, < fj < 1 is the occupation fac- 
tor of the j-th site, and Dkj = Ek — Ej — F(xk — xj). The 
net current from j to k is given by Zjy = efTjy — r*). The 
transport problem amounts to solving the current conserva- 
tion equations ^k^jhj = for ffs at a given fixed F. One 
conventional way to present the result is in the form fj = 
{exp[(e 7 - — fjj)/ksT] + where /jj is termed the chemi- 

cal potential of the j-th site. In the most relevant situation, 
A kj = \Ej-fii\ + \Ek - fJk\ + \Ej - E k \ > k B T, the currents I kj 
are given by 1 14] 

( 2 1 i A kj \ ■ u lj-lk n . 

Ikj - «fi*p[--\*-xj\ ~ ) sinh ^-' P) 

where Z,j = fij — F x i i s the electrochemical potential. Know- 
ing /j ; 's (or //s), one then computes the total current I and the 
resistivity p(F,T) = el/F. 

Ohmic regime. — In general, the above transport equations 
are nonlinear in /jj but in the Ohmic limit, F — » 0, T > 0, 
both fjj and £,j — ^ are smaller than T and Eq. (0 becomes 
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hjRkj, 



where — Roexp[(2/a)\xk 



■Xj 



(Ay/UcsT)]. Therefore, the system is equivalent to the net- 
work of resistors R^, while t, ,-/e is simply the voltage at the 
j-node OJ01. Here /?o = kgT / (e 2 y) and A» is to be computed 
at fjj = [ifr = 0. The resistances differ by many orders of 
magnitude from one another and this is what makes the prob- 
lem nontrivial. Let us review the basic ideas developed for its 
solution. A crucial step 111 |2|] is to realize that among many 
possible current paths through the network, a single optimal 
one heavily dominates the conduction. In ID this path has the 
topology of a linear chain, so that we can label the nodes of 
the path consequtively: 1,2, ...,/,/+ 1, etc. The typical sepa- 
ration in position (xm) and in energy (Em) between these nodes 
follows from the conventional Mott's argument |l|,|2j]: 




£m(T) _ x M (T) 
kr>T a 



■ um : 



, k B T ; 



1 

ga' 



(4) 



Thus, the logarithm of the resistance of a typical link is 
ln^j+i ~ xm /a + ZM/ksT ~ um- In addition to typical links, 
there are places in the conductor where hopping sites are lo- 
cally absent in a large area A of the (x,e)-plane. The probabil- 
ity of every such "break" is very small, ~ exp(-gA), but its 
resistance far exceeds that of a typical link. A careful analysis 
of Raikh and Ruzin (RR) [5] showed that the dominant contri- 
bution to p is provided by breaks that are shaped as diamonds 
and have the sizes u&a/2 and 27b in the position and energy 
dimensions, respectively, see Fig. da) LJ. Here u a = T /T. 
These "optimal" breaks have resistances R& — Rq exp(i<n) and 
are present in concentration Nq ~ exp(— uq/2)/xm, which 
yields p - Rq,Nq, i.e., |5] 



FIG. 1: (a) Ohmic and (b) non-Ohmic breaks. The dots represent 
hopping sites. The horizontal thin lines are the links of the optimal 
path carrying the current across the breaks. The width of each break 
is auj/2. 



the electrochemical potential difference across the z'-th link of 
the optimal path. Since xm is the correct estimate of the typi- 
cal length of the link, we have F ~ (A^,-) /xm, where the angle 
brackets denote averaging over i. Introducing P^, the proba- 
bility distribution function (PDF) of the random quantity A^i, 
we can write this as 



F~— /^(A^)A^. 

xm 



(7) 



lnp(0,r) ~T /2T, r<7b. 



(5) 



Let us establish the functional form of P^. The behavior of 
P^(At,) at At, <^ksT is determined by links with Ohmic re- 
sistances Roe 11 , where u <C u\. They remain operating in the 
Ohmic regime, and so for them A^j/fcjjT = exp(« — ui) -C 1. 
Using RR's results for the PDF of the (diamond-shaped) 
Ohmic breaks, we obtain 



As customary in this field, we will not attempt to compute the 
pre-exponential factor of the resistivity 111 211 . 

If the ID Ohmic VRH seems like a nontrivial problem, it 
becomes even more so in the non-Ohmic regime we wish to 
discuss next. In particular, we address the following new set 
of questions: What is the highest F when Eq. still applies? 
Do breaks continue to play a role at larger Fl How does p 
depend on F at such fields? Our answers to these questions 
are given below. 

Intermediate fields. — We begin by showing that the region of 
validity of Eq. (0 is very limited. Indeed, in the derivation it 
was assumed that the voltage drop across each link is less than 
ksT/e. However the voltage drop on every optimal break is 
equal to IR& = (F /ep)/fo —F/eNa; therefore, Eq. (jSJi is jus- 
tified only at F <C Fq ~ (ksT /xm)sxp(—Tq/2T). To compute 
p at larger F a different approach is needed. Our strategy re- 
sembles that of Ref . 1 1 3 ] . We assume that a given fixed current 
I flows through the system and compute F = (d%/dx). Then 
we find p from the relation 



C T 
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At this stage we focus on the intermediate field (current) 
regime, defined by «q. Let A^,- = — be 



where 1 <C \n(ksT / At,) < uj — um and C is a prefactor with 
a subleading dependence (presumably, a logarithmic one) on 
At, and T. Note that the typical links, u ~ um, are included in 
Eq. © [they correspond to the lowest At, ~ k B T exp(uM — «/) 
at which this formula is still valid]. Consider now the links 
with At,i ^> k B T, which are necessarily atypical, i.e., some 
sort of breaks. The geometry of such non-Ohmic breaks can 
be established as follows. In a vicinity of a non-Ohmic break 
the largest voltage variation is across this break; therefore, the 
electrochemical potential has nearly the same value, say, t,[, 
on the sites adjacent to the left boundary of the break and an- 
other value, t,j + i for those on the right side. Assuming that 
At,i ~ A/j,, (borne out by the results) we see from Eq. (0 
that the current through a pair of sites j and k located on 
the opposite sides of the break is Iy ~ eyexp[— (2/a)\x/ c — 
Xj \ — (Aiti/lksT) + (Afii/ksT)]. The shape of the non-Ohmic 
ui = In ( —\ . (6) break should be such that the inequality Itj < / is obeyed. The 
optimal, i.e., the most probable breaks are those of minimal 
area. Using only elementary geometrical considerations one 
readily finds that the optimal shapes comprise a continuum 
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of curvilinear hexagons of area A; = au](k B Tui + A^,)/2, see 
Fig. Hlb). Hence, /^(A^) ~ exp(-gA,) for > k B T, or, 
more accurately, 

C ( Tu 2 , m/AE \ 

Equations (jSJl and (|9jl match by the order of magnitude at St, ~ 
k B T. Substituting them into Eq. 0, we find the the main 
voltage generators in the system are the non-Ohmic breaks 
with Acj ~ ksTo/ui. For F we obtain 

F ~ (Ck B T Q /auj)exp(-uj /2u 2 M ). (10) 

Solving this for uj and using Eq. (|6j, we arrive at 

llT Q fvk B T\ k B T J, k B T 

where v = Cuj/u^. The main message of this analysis is that 
atypically large resistors are progressively eliminated as F in- 
creases (uj decreases). Note that at uj ~ uu the distinction 
between breaks and typical links disappears, so that the result 
Acj; ~ k B To/ um for the breaks is a good estimate of the overall 
(At,). Thus, at this point F ~ k B T /a and therefore v ~ C ~ 1. 
Strong fields. — A quick if unilluminating derivation of p at 
F ^> k B T /a assumes that Eq. (1101 holds throughout the non- 
Ohmic regime. In strong fields, where uj <C «m, the exponen- 
tial in Eq. (II 01 drops out leaving one with F ~ Ck B To/auj. 
Combined with Eq. (|6}, this entails Eq. (0. Notably, this 
formula differs from Eq. (4.29) of Ref. [11], by a large log- 
arithmic factor in the exponential. In view of this discrep- 
ancy (which has a conceptual importantce, cf. the introduc- 
tion), we present a rigorous derivation of Eq. Q that clarifies 
the physics involved and yields the exact result C — 8 for the 
coefficient C. 

In strong fields only the forward hops need to be consid- 
ered J^Eal, so that for the r'-th link of the optimal path we 
have / = eTf+n or, equivalently [cf. Eq. Q] 

fi{\ - /i+i)exp[K/ - (2/a)\xi+\ = 1, (12) 
e,>e i+ i, EiEEEi-Fxi. (13) 

A necessary condition for the existence of a physically accept- 
able solution, < /, < 1, of Eq. H2\ is 

< x i = oui/2 (14) 

(thus, the typical hopping length is field-dependent). A suf- 
ficient condition is known from the theory of a Disordered 
Asymmetric Simple Exclusion Process fl7ll . which is de- 
scribed by the identical equations. It is just a bit more strin- 
gent, \xi+\ —Xi\ < xj — aln2. Since we are interested in the 
case ui ^> 1, it is a legitimate approximation to assume that 
the condition dl4> is both necessary and sufficient. Only the 
existence of the solution matters, not the actual values of /,-. 
(Those will be commented on at the end.) We expect on phys- 
ical grounds that the optimal path does not deviate indefinitely 




FIG. 2: Plots of E(x) (thin line), E < (x) (thick line) and the optimal 
path (dotted line). Ticks and dots represent lattice sites and their 
energy levels, respectively; xj = 9b is assumed. 

from the equilibrium chemical potential, sup |e,-| < °°. Then 
the standard percolation theory argument II 1411 entails 

F = minlim[(ei —£;)/(*,•— *i)]. (15) 

Indeed, if F were smaller than the right-hand side of this equa- 
tion, then the contiguous path that obeys the constraints ( II 31 
and (1141 would not exist. If F were larger, the chosen path 
would be short-circuited by another one. To compute F from 
Eq. (1151 we use the fact that statistics of the random points 
(xj,£j) is again an uncorrected Poisson distribution of den- 
sity g. 

Consider the following lattice version of the problem. Let 
all coordinates be restricted to integer multiples of b <C x\. On 
every site xj — jb let there be a discrete set of energy levels 
E n (xj), n € Z, drawn from an independent Poisson distribu- 
tion with the average energy spacing 1 jgb. Suppose next that 
{x\ , Eq (x\ )} is the starting (the leftmost) site of the sought op- 
timal path. Denote by E(x) the largest E„(x) reachable from 
this site by a path that obeys the constraints ( 1131 and (1141 . then 
F = lim x ^„o{[eoC>i) -E(x)]/(x-xi)} = -1im x ^ 00 [E(x)/x]. 
Function E(x) satisfies the obvious recurrence relation: E(x) 
is equal to the largest £ n (x) that still does not exceed E K (x), 
where E < (x) = maxE(y), x — xi < y < x. This relation is il- 
lustrated graphically in Fig. |2] It is easy to see that function 
£<(x) is step-like, i.e., that there is an increasing sequence 
of coordinates Xi, such that E<(x) = E{xt) = const for xi < 
x < Xi+\. Figure [2] aids in observation that given xi, the next 
site in this sequence, Jtj+i, is simply the one with the largest 
E(y) less than £<(x,) among all the points that belong to the 
vertical strip Xj < y < Xj +xj. These x, are in fact the coordi- 
nates of the optimal path (more precisely, they become such 
once one further optimizes with respect to the path's starting 
point). By virtue of this property, AEj = E < (xi) — E < (xi+ 1) 
and Axj = xi+\ — x, are uncorrected random numbers with av- 
erage values 2k B To/ui and aui/4, respectively (in the limit 
b <C xj). This entails F = (AE,)/(Ax,) = Ck B To/auj with 
C = 8. 

It is worthwhile to clarify why Eq. is not invalidated by 
inevitable large voids in the distribution of sites in the (x,e)- 
plane. Whenever the optimal path meets such a void, A£, be- 
comes much larger than its average value. However, the corre- 
sponding contribution to F grows linearly with the size of the 
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void whereas the probability of the void drops exponentially. 
These rare large voids have thus a negligible effect on F(uj). 

The derivation above indicates that the site occupations in 
the energy band |e,-| < ksTo/uj ~ (kgToFa) 1 ' 2 that carries 
the current, deviate so much from their equilibrium values, 
ln(l jfj — 1) = Ej/IcbT, that T plays no vital role in the trans- 
port. This agrees with the early result of Shklovskii 11611 
that at large |e,-| the real temperature T in the last equation 
gets replaced by the "effective" one, T e g = Fa /2kg. How- 
ever, at e/s of relevance to VRH there exist large correc- 
tions, fcgr e ffln(l//) — 1) — £j ~ {ksToFayl 2 , to the asymp- 
totic large-e behavior that prohibit one to derive Eq. Q via 
a simple replacement T — > T e ff in the weak-field formula, 
Eq. 

Experimental implications. — The length L of real ID con- 
ductors is seldom very large. Because of that their resistiv- 
ity in a VRH regime fluctuates strongly from one sample to 
the next, depending on what configuration of breaks is real- 
ized 1 4, S\. Here we discuss only the ensemble average value 
of In p. For the Ohmic case, in the most common situation 
1<CK Tq/T, where L = \n{L/xM), the result for this quan- 
tity is known to be H (lnp(0,T)) ~ y / 2LT /T. In this 
regime practically all the applied voltage V = FL/e drops on 
the single largest break |5|. Hence, the Ohmic behavior lasts 
as long as eV < k B T . In some range of eV above k B T, the 
same break remains the bottleneck for the transport. Using 
Eq. (0 and averaging over the position of the break along the 
chain, we obtain the exponential dependence, 

ln(p/p(0,r)) ~ -eV/2k B T, k B T < eV <eV L . (16) 

The upper limit Vi where this formula is valid can be esti- 
mated from the condition eVr, ~ (At,) ~ k B To/uj, which gives 
eV L ~ k B (TT / ' L) x l 2 . At V > V L the infinite-chain results, 
Eqs. (till and Q, apply. 

On a qualitative level, existing experimental data are con- 
sistent with our theory. For example, a characteristic S-shaped 
log/ — logV curve KB follows naturally from our equations. 
To unabiguously verify them an earnest effort and a thoughtful 
experimental design are required. At present, our knowledge 
and degree of control over disorder and other parameters of 
ID systems remain poor; therefore, a wide dynamical range 
of data is needed to compensate for many uncertainties. The 
most suitable for this purpose are systems where localization 
length is not too short, e.g., high-quality carbon nanotubes. 
Their resistance has to be measured over a broad range of low 
T, where p is determined by the exponential terms we com- 
puted rather than some power-law prefactors we left undeter- 
mined here. A care should be taken to remain in the validity 



domain of the theory, L S> xm ~ a y/To/T. The parasitic ef- 
fects of contact resistance la |8|] need to be done away with 
reliably. If the system is quasi- ID rathen than ID, hopping 
across the chains must be negligible. Finally, a careful 
statistical analysis of the data is likely to be needed I HI. 

An important qualitative prediction of our theory is that 
enormous sample-to-sample variability of the low-T hopping 
resistance becomes less of an issue in the non-Ohmic regime. 
This may facilitate both the experimental studies we call for 
and the future theoretical work that would be needed to in- 
clude the effects of Coulomb interaction into the model. Such 
coupled experimental and theoretical studies will likely to 
shed much light on the intriguing transport properties of low- 
dimensional conductors. 
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